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A local formulation of lattice QCD without unphysical fermion zero 
modesB 
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The problem of unphysical zero modes in lattice QCD with Wilson fermions can be solved in a clean way by 
including a mass term proportional to iip-y^T^tp in the standard lattice theory with Ni = 2 mass degenerate Wilson 
quarks. We argue that up to cutoff effects, this lattice theory is equivalent to standard lattice QCD, for suitable 
choices of the mass parameters and with a natural re-interpretation of observables. On-shell 0(a) improvement 
can be implemented in a straightforward way. 
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1. INTRODUCTION 

Lattice QCD with Wilson quarks provides an 
attractive framework for non-perturbative com- 
putations of many interesting phenomcnological 
quantities. Its main disadvantage consists in the 
explicit breaking of all axial symmetries. As 
a consequence the Wilson-Dirac operator is not 
protected against zero modes even at finite values 
of the quark mass. While this is not a theoreti- 
cal problem for the full theory, it severely limits 
the applicability of the quenched approximation 
as the occurrence of unusually small eigenvalues 
may lead to large fluctuations in observable quan- 
tities. In practice this means that one cannot di- 
rectly treat the light quark masses but has to ex- 
trapolate results obtained with relatively heavy 
quarks. It is important to appreciate that this 
limitation is completely independent of the usual 
requirement of not too light pion masses, dictated 
by the finite space-time volume used for the sim- 
ulation. Indeed quite often quenched simulations 
are limited by the occurrence of "exceptional con- 
figurations" at quark mass values where finite vol- 
ume effects are still completely negligible. 

Various recipes have been proposed to deal 
with this problem [1-3] but to our knowledge no 
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completely satisfactory solution has appeared in 
the literature. In order to fill this gap, we consider 
the lattice action first introduced by Aoki |Q. 
More precisely, if Dw denotes the (0(a) im- 
proved) Wilson-Dirac operator for a quark with 
bare mass mo, we consider the fermionic action 
for Nf — 2 quarks with the lattice Dirac operator 



-Dtwist — 



(1) 



The additional term has a non-trivial flavour 
structure (the Pauli matrix r 3 acts in flavour 
space), and will be referred to as a chirally twisted 
mass term. It is obvious that this lattice Dirac 
operator is protected against zero modes for any 
finite value of fj, q , 



det Awist = det(D^ w D w + fi 2 ) > 



(2) 



To use this lattice action as a regularization of 
"exceptional configurations" has therefore been 
proposed in ref. |J, implying, however, that the 
mass parameter /x q be extrapolated to zero at the 
end. As will become clear below, this is not nec- 
essary. It can be shown that lattice QCD with a 
chirally twisted mass term (twisted lattice QCD) 
is, up to cutoff effects, completely equivalent to 
standard lattice QCD, and thus provides an al- 
ternative regularization of the massive theory. 

This article is organised as follows: we first take 
a closer look at the classical continuum limit of 
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twisted lattice QCD (sect. 2). Assuming univer- 
sality we sketch an elementary proof that twisted 
and standard lattice QCD are equivalent in the 
continuum limit (sect. 3). In sect. 4 we determine 
the possible counterterms, and briefly introduce 
our choice of renormalization scheme. We then 
discuss 0(a) improvement (sect. 5) and we end 
with a few comments on further applications. 

2. CLASSICAL CONTINUUM THEORY 

To get acquainted with twisted QCD it is in- 
structive to have a closer look at the correspond- 
ing classical continuum theory, described by the 
lagrangian 

C f (x) =lp(x) (p+m + ifiq^T 3 ) ip(x). (3) 
Defining the usual isospin currents and densities, 



(4) 
(5) 
(6) 



the partial conservation of axial and vector cur- 
rents reads □ 



d^Al = 2mP a , 



1,2 



dpAl = 2mP 3 + in q iJi{j, 

d^V? = -2^e ab P\ a, b -1,2 



(7) 
(8) 
(9) 
(10) 



A global chiral (non-singlet) rotation of the fields, 

= exp(ia75r 3 /2)'0, 

^' = i>exp(iaj 5 T 3 /2), (11) 

with tana = fi q /m q transforms the Lagrangian 
to its standard form, 

(12) 



C' f {x)=^\x) (p+m')i;'(x), 

with the quark mass m! = 
isospin currents and densities in the primed basis 



i 2 + nl 



The 



2 We define the antisymmetric tensor e a b with indices 
a, b = 1,2 and normalization ei2 = 1. Summation over 
repeated indices is always understood. 



are related to the original fields through (a, b 
1,2) 



A'l = cos(a)A° + sin(a)e ab V;, 



A . 



r ,a 



A 3 



V'l = cos(a)V; a + sin(a)e afc ^, 

V' 3 = V 3 

pia _ pa 

P' 3 = cos(a)P 3 -Hsin(a)|M 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



and satisfy the PCAC and PCVC relations in 
their standard forms, i.e. 

d^A'l = 2m'P'\ dpV'l = 0, a =1,2, 3. (19) 

Similar relationships can be worked out between 
more complicated composite fields , and we thus 
conclude that the classical twisted theory can be 
mapped back to the standard formulation by a 
simple axial rotation of the fields, with a unique 
rotation angle which is fixed by the ratio of the 
mass parameters. Both formulations are thus 
equivalent and contain the same physical infor- 
mation. 

3. BEYOND THE CLASSICAL THEORY 

In order to establish the equivalence between 
twisted QCD and the standard formulation of 
QCD beyond the classical continuum limit, it is 
convenient to make use of a regularization which 
does not break the flavour symmetries. This is 
possible Hf?} if the lattice Dirac operator satisfies 
the Ginsparg- Wilson relation S , 



£>75 + 75-D = aD^D. 



(20) 



An explicit solution for D has been given by Neu- 
berger but in the following we only need to 
know that D is a local operator which satisfies 
the Ginsparg- Wilson relation and has the conju- 
gation property D 1 ' — 75D75. It then follows that 
the matrix [p~0| , 



75 = f 75(1 ~ aD), 



(21) 



is hermitian and unitary. Therefore it may be 
used to define left handed quark fields through 



3(1 -w, 



^=^(1 + 7*), (22) 
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and right handed fields are obtained with the 
complementary projectors. Using these chiral 
fields a lattice regularized version of twisted QCD 
is specified by the action 



Sf = 



[^lD^l + V'r^r 



(23) 



In the chiral limit m = /i q = this action has 
an exact SU(2) x SU(2) flavour symmetry, and in 
particular the axial transformation [cf. eq. (fTl])], 

ip' L = exp(-zar 3 /2)i/'L, 

ij/ h = ip L exp(iaT 3 /2), 

ip R = exp(iar 3 /2)^ R , 

^R = V' R exp(-iar 3 /2), 



(24) 



is part of the symmetry group. Therefore, this 
transformation leaves the form of the lattice ac- 
tion invariant and merely transforms the param- 
eters, 



m = m cos(a) + /i q sin(a), (25) 
/i q = — m sin(a) + /i q cos(a), (26) 

where the choice tana = /i q /m again produces 
the standard action with /i q = and m' = 
(to 2 + /i 2 ) 1 / 2 . Since the flavour symmetry in 
the basis of chiral fields takes exactly the same 
form as in the continuum, we may also reproduce 
eqs. (13-19). More generally, bare operators in 
a given representation of the flavour symmetry 
group are rotated covariantly by the same an- 
gle a. Finally we note that the measure of the 
functional integral is left invariant under the axial 
(non-singlet) transformation, eq. (|24|). We have 
thus found a regularization where (at least in a fi- 
nite space-time volume) the bare parameters and 
correlation functions of twisted QCD are related 
to the ones of standard QCD through a change 
of variables in a well-defined functional integral. 
This yields exact identities at any fixed value of 
the lattice spacing, and gives a precise meaning 
to the equivalence between the two theories. For 
example, if O k and 0' k , are, for k = 1, . . . , n, bare 
gauge invariant composite operators, which are 



related through the axial transformation (|24|), we 
obtain the identity between the (bare) connected 
correlation functions, 



(Oi(xi) 
(0[( Xl ) 



D ( t "\\ twisted QCD 

' 0„ {x n )) c, 



(27) 



provided the bare parameters /i q , to and m! are 
related as indicated previously. 

In order for these relations to carry over to the 
continuum limit, it is convenient to renormalize 
both mass parameters, /i q and m in the same way, 
such that the ratio /iq/m, and thus the angle a, 
remains unchanged. If the chosen renormaliza- 
tion scheme respects the global flavour symme- 
tries, then the relations between correlation func- 
tions of composite operators are exactly the same 
as in the bare theory, and thus also hold in the 
continuum limit. Mass independent renormaliza- 
tion schemes obviously have this property and 
in addition preserve relations between renormal- 
ized parameters, but it turns out that a much 
wider class of renormalization schemes is admis- 
sible 0. 

In lattice QCD with Wilson fermions an ax- 
ial transformation does not correspond to an ex- 
act symmetry. However, assuming universality 
of the continuum limit, we may conclude that 
the twisted and standard formulations of lattice 
QCD with Wilson fermions are indeed equivalent 
up to cutoff effects. In fact, this can be proven 
rigorously to all orders in the loop expansion, by 
exploiting reparametrization invariance and the 
Quantum Action Principle p]| . In particular, 
this proof holds for any regularization for which 
a power counting theorem can be established, re- 
gardless of whether or not the flavour symmetries 
are broken by the regularization. 

4. RENORMALIZATION 

On the lattice with Wilson fermions only the 
isospin SU(2) symmetry remains exact, and the 
flavour symmetry is further reduced by the in- 
troduction of the twisted mass term. This is 
in contrast with the lattice regularization based 
on Ginsparg- Wilson fermions, where the symme- 
tries are exactly the same in both (standard and 
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twisted) formulations, albeit realized in a differ- 
ent way. 

A careful analysis shows that, apart from gauge 
invariance and invariance under discrete four- 
dimensional rotations, the following symmetries 
remain exact on the lattice: 

• a global U(l) symmetry corresponding to 
fermion number conservation, 

• a residual U(l) flavour symmetry with gen- 
erator T 3 /2, 

• charge conjugation, 

• Pp: parity combined with flavour exchange, 

• F: flavour exchange combined with /i q 

Furthermore, it is possible to show that reflection 
positivity |l2) still holds. 

Taking into account these symmetries, twisted 
lattice QCD is renormalizable by power counting. 
The only new counterterm of dimension < 4 is 
proportional to ifiqip^T ip, which leads to a mul- 
tiplicative renormalization of fi q . In particular we 
note that a term of the form e^ upa tr {F^ v F pa } 
is excluded by the symmetry Pp, whereas in- 
variance under F excludes an additive renormal- 
ization of /i q (i.e. a counterterm proportional to 

We would like to impose renormalization con- 
ditions such that the relation between renormal- 
ized lattice QCD and its twisted version is as sim- 
ple as in the classical continuum limit. However, 
due to the reduced flavour symmetry with Wilson 
fcrmions some preparation is needed. In particu- 
lar, one has to determine the relative normaliza- 
tion between composite operators that belong to 
the same representation of the continuum flavour 
symmetry, but are not related by lattice symme- 
tries. This can be achieved by imposing contin- 
uum chiral Ward identities as normalization con- 
ditions for the lattice operators (ll| . An example 
is provided by the isospin currents: while the con- 
served vector current is protected against renor- 
malization, the correctly normalized axial current 
(A R )^ has to be determined by imposing current 
algebra relations as normalization conditions JTij . 



We then expect that relations like eqs.( 13-18) 
hold for the correctly normalized currents, pro- 
vided the angle a is defined using the Ward iden- 
tity masses. More precisely, we define a bare cur- 
rent quark mass m through (some suitable matrix 
element of) the PCAC relation, 

d^(A n )l = 2mP\ (28) 

and define a through 

tana — fi q /m. (29) 

Note that m differs from the subtracted bare 
quark mass m q = mo — m c . In the renormalized 
theory, we thus fix the ratio of the mass renor- 
malization constants 

m R = Z m m q , fi R = Zpfiq, (30) 

such that it coincides with the slope at the origin 
of the function m(m q ). 

5. O(a) IMPROVEMENT 

Given the symmetries of the lattice theory it is 
a straightforward exercise to list the possible 0(a) 
counterterms to the action. The procedure has 
been explained in great detail in ref. ||l5f and ap- 
plies to the case at hand, taking into account the 
modified set of symmetries given in the previous 
section. We expect the renormalization and 0(a) 
improvement of the massless theory to be the 
same as in standard lattice QCD [13-15]. This 
follows from analyticity of correlation functions 
in the mass parameters 77i q and /i q , which is how- 
ever guaranteed only with an infrared cutoff, e.g. 
in a finite space-time volume. 

With this assumption, on-shell 0(a) improve- 
ment of the action is achieved by using the 
Shcikholcslami-Wohlcrt improved Wilson-Dirac 
operator Jl6|] , and the following parameterization 
of the 0(a) improved renormalized coupling and 



quark masses, 

9k = .9o z g(SoW), (31) 

m R = m q Z m (gl,ay) 7 (32) 

MR = Mq-Mi?0 ( 33 ) 

.9o = 9 2 (l + b g am q ), (34) 

m q = m q + 6 m am q + o m a/i q , (35) 
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H q = Mq(l + b^arriq), 
m q = mo — m c . 



(36) 
(37) 



This notation already implies our choice of a 
mass-independent renormalization scheme [ |l8[ 
p^[ . 0(a) improvement of the action thus intro- 
duces the improvement coefficients & M and b m , in 
addition to the standard coefficients b m and b g . 

As for the improvement of operators we here 
give a few examples which appear in the PCAC 
and PCVC relations. Following ref. we define 
the improved operators (a = 1,2), 

(Aj)^ = Al+c A adf,P a + ati~bAe ab V*, (38) 

{ViYu = V" + c v ad v T«+ati q b v e ab A b fl , (39) 



(Pi) a = P a , 



(40) 



and obtain the renormalized improved operators 
through (a = 1,2): 



2 A (l + i) A am q )(A I );, 

Z v (l+6 v am q )(K)^ 
Zp(l + fe P am q )(F I ) a . 



(41) 
(42) 
(43) 



Note that we may avoid technical problems asso- 
ciated with power divergences by restricting at- 
tention to the first two isospin components of 
these operators. This can be done without any 
loss of information due to the expected restora- 
tion of all chiral Ward identities in the continuum 
limit. 

With these definitions the 0(a) improved angle 
a is defined by 



tana = 



Mr. 



Z m rh a 



(44) 



where /x R and m R are obtained from matrix ele- 
ments of the renormalized 0(a) improved PCAC 
and PCVC relations. Again this requires that the 
ratio Z fl /Z m be fixed as explained in the previous 
section. 

In order to check the general framework of 
renormalization and 0(a) improvement described 
above, we carried out a perturbative calculation 
using on-shell correlation functions derived from 
the Schrodinger functional |2(J . The details of the 
calculations are very similar to the studies of the 
standard theory (2|^^] and will be described else- 
where [frii . As a result we confirm the theoretical 



expectations. In particular, all improvement coef- 
ficients are found to be functions of the bare cou- 
pling alone, and do not depend e.g. on the angle 



a. To lowest order we find b 



(0) 



r(o) 



,(o) 







and b 



(0) 



-i, and the corresponding one- loop 
coefficients are available, too [ [rT| . 

6. CONCLUSIONS AND OUTLOOK 

There remains little doubt that twisted lattice 
QCD and standard lattice QCD for Nf = 2 mass 
degenerate quarks are equivalent theories in the 
continuum limit, up to 0(a 2 ) effects if 0(a) im- 
provement is implemented. Since twisted lattice 
QCD is protected against unphysical fermion zero 
modes, the problem of exceptional configurations 
has thus found a clean solution. Furthermore, 
realistic studies of QCD with Nf > 2 can be eas- 
ily done by using the (improved) standard lattice 
action for the heavier quarks. 

In full QCD the problem of exceptional con- 
figurations is absent in principle. However, it 
has been argued that close to the chiral limit, 
standard algorithms may still experience techni- 
cal problems, showing up e.g. as a poor sampling 
of configuration space p3|]2~i| ] . Twisted QCD may 
therefore be useful in this context, too. Inciden- 
tally, full QCD simulations using (unimproved) 
twisted lattice QCD have already been carried 
out, albeit with a different motivation pq |. 

In view of 0(a) improvement, there are pre- 
ferred choices of a. If a is small, perturbative 
estimates of the new &-coefficients are sufficient, 
while large fluctuations due to close zero modes 
are already excluded. Another attractive possi- 
bility is the choice a — tt/2, such that the physi- 
cal quark mass is entirely defined in terms of /i R . 
In this case all the standard 6-cocfficients are not 
needed, and, in particular, c/q = go- This latter 
fact may be interesting for full QCD simulations, 
as continuum and chiral extrapolations may be- 
come easier. 

Since theories with different values of a are not 
related by a lattice symmetry, operator mixing 
problems may look quite different. A simple ex- 
ample is provided by the iso-singlet scalar density 
where the usual additive renormalization at a = 
can be avoided by working at a = tt/2. 
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Dynamical fermion simulations being very ex- 
pensive it has been proposed to obtain informa- 
tion about full QCD by extrapolating from the 
unphysical region where Nf < [26-28]. For 
this "bermion" approach the problem of close zero 
modes is much more severe than in the quenched 
approximation, and may again be solved using 
twisted lattice QCD. 

Originally, Aoki introduced the chirally twisted 
mass term to discuss the phase structure of the 
standard lattice theory Q . At strong coupling he 
established the existence of a phase with spon- 
taneous breaking of parity. We do not exclude 
that such an Aoki phase also exists at weaker 
couplings. However, even if it does, this would 
be relevant only in large space-time volume and 
very close to the chiral limit (mass parameters of 
0(a 2 )) where simulations are not practical any- 
way (cf. ref. ||). 

In conclusion, we believe that twisted lattice 
QCD provides a very promising alternative for- 
mulation of lattice QCD. 
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